Optical conductivity of rattling phonons in type-I clathrate Ba 8 Gai 6 Ge 3 o 
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A series of infrared-active optical phonons have been detected in type-I clathrate BagGai6Ge3o 
by terahertz time-domain spectroscopy. The conductivity spectra with the lowest-lying peaks at 
1.15 and 1.80 THz are identified with so-called rattling phonons, i.e., optical modes of the guest ion 
Ba 2+ (2) with T\ u symmetry in the oversized tetrakaidecahedral cage. The temperature dependence 
of the spectra from these modes are totally consistent with calculations based on a one-dimensional 
anharmonic potential model that, with decreasing temperature, the shape becomes asymmetrically 
sharp associated with a softening for the weight to shift to lower frequency. These temperature de- 
pendences are determined, without any interaction effects, by the Bose-factor for optical excitations 
of anharmonic phonons with the nonequally spaced energy levels. 

PACS numbers: 63.20.Ry, 78.30.-j, 82.75.-z 



During the past decade, thermoelectric materials such 
as clathrates and filled-skutterudites have renewed an in- 
terest in phonons These conductors are formed by 
polyhedral building blocks, where each polyhedral cage 
accommodates a guest ion like an alkali-earth or rare- 
earth element. If the cage is oversized, the guest ions 
can vibrate with large amplitude around an on-center 
or off-center site in the cage potential. These Einstein- 
like local modes have been called rattling phonons. The 
renewed interest above has been paid to interactions of 
these rattling phonons with acoustic phonons propagat- 
ing through the cage network and carrying heat entropy, 
and, equally or more interestingly, with conduction elec- 
trons. However, these issues are still far from being well 
understood that even the charge dynamics of a single ion 
in rattling vibration have been little known. 

This paper reports on the optical conductivity of rat- 
tling phonons detected in a type-I clathrate BasGai6Ge3o 
(hereafter abbreviated as BGG), featuring the anhar- 
monicity effects on the conductivity spectra in compari- 
son with model calculations. 

This compound belongs to a family of ternary type- 
I clathrates having cubic crystal structure with space 
group The unit cell of the host framework 

consists of 6 tetrakaidecahedra and 2 dodecahedra. The 
latter smaller cages around 2a sites occupy the body- 
centered- cubic sites, while the former oversized cages 
around 6d sites line up along the principal axes by sharing 
both the pentagonal faces with the dodecahedral cages 
and the hexagonal faces with neighboring tetrakaideca- 
hedra. Every cage encapsulates one Ba 2+ ion, and the 
guest ions in the smaller (2a) and oversized (6d) cages 



are defined as Ba 2+ (1) and (2), respectively. These guest 
ions satisfy the Zintl rule for charge compensations that 
the more electro-positive guest atom donates electrons 
to the more electro-negative cage; 16 electrons are trans- 
ferred from 8 Ba atoms to 16 Ga atoms of the cage. With 
finely tuned Ga/Ge concentration, therefore, the system 
can be a heavily doped semiconductor, in general, having 
both the charge-carrier's sign and density controlled. 

So far the electrical resistivity, thermal conductivity 
and specific heat measurements have been made on BGG 
having different carrier densities of n- or p-type carriers 0, 
4]. The nuclear density map for the Ba 2+ (2) ion in the 
oversized cage obtained by neutron diffraction Q indi- 
cates that Ba 2+ (2) resides on center, while other ions, 
Sr and Eu, move away to off-center sites, being sup- 
ported by the band structure calculations^. On the 
contrary, recent more detailed structural analyses Q indi- 
cate that the Ba 2+ (2) ions are slightly displaced from the 
center depending on temperature and also carrier-type. 
This structural indication seems to be consistent with 
Raman scattering measurements of low lying additional 
modes 0- 

Single crystals of p-type BGG are grown by a self-flux 
method Q. To obtain transmitting signals as strong as 
possible through the conductive sample, single crystal 
disks of 2- 5mm in diameter and 0.5 -1.0 mm in thick- 
ness are polished down to about 20 y«m in thickness 
using diamond abrasive films. Terahertz time-domain 
spectroscopy (THz-TDS) measurements covering the fre- 
quency range of 0.2 - 3.5 THz (0.8 - 14.5 meV) are carried 
out with the spectrometer (RT-20000, Tochigi Nikon Co. 
Ltd) which uses a standard technique for the transmis- 
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FIG. 1: (color online) (a) Time-evolution of THz-wave elec- 
tric fields obtained from three different sets of measurements; 
vacuum, the substrate only and the type-I BasGaieGeso 
(BGG) sample glued onto the substrate and (b) the Fourier- 
transformed power spectra |i?(u;)| 2 . 
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FIG. 2: (color online) Temperature-dependent conductivity 
spectra a(uS) = cri(u)) + i<72(w) of BGG. For convenience, 
ITHz = 33 cm' 1 = 4.1 meV = 48 K. 



sion configuration, see, e.g., our preceding paper[9[. 

Figure [TJa) shows the time-evolution of transmitting 
electric fields through the vacuum, the substrate only 
and the sample glued onto the substrate, while Fig. QJb) 
shows the Fourier-transformed power spectra. Then the 
refractive index spectra are determined by taking into 
account multiple reflections at the surface and boundary 
between the sample and substrate 

Figure [5] shows the conductivity spectra <j{oj) = 
o\ (u>) + icr-i (u>) . The real part cti (u>) takes several distinct 
peaks, while the imaginary part 02(0;) shows correspond- 
ingly the derivative with frequency. These discrete spec- 
tra, reminiscent of Lorentz spectra for optical phonons, 
are superimposed on the monotonous background con- 
ductivity. The magnitude of the background at room 
temperature is about 10f2 _1 cm _1 comparable to the dc 
conductivity and, with decreasing temperature, decreases 
almost down to zero, again consistent with the temper- 
ature dependence of the dc conductivity [8j. Therefore 
it is concluded that the conductivity spectra consist of 
both an almost constant Drude-like contribution due to 
the semiconducting carriers and the Lorentz-like spectra 
from optical phonon modes. 

For the assignment of thephonon modes, we have made 
first-principle calculations [7| based on the density func- 
tional approximation for infrared-active optical modes of 
T\ u symmetry. Table U lists the phonon frequencies la- 



TABLE I: Calculated phonon frequencies and fitting parame- 
ters (6.6 K) of BGG. Due to the broad spectra around 2.4 THz 
(Fig. [2]), the analysis for v% and Va is made with a single 
Lorentz curve (Eq.[2}. 
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0.96 


1.15 
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8.8 


v% 


Ba(2) 


1.85 


1.80 
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2.70 
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0.4(1) 




cage 


3.05 


3.30 


0.07 


0.6(4) 



beled as vi(i = 1 — 6) together with fitting parameters 
at T = 6.6 K. The low- lying peaks observed at 1.15 and 
1.80 THz can be assigned as the modes labeled v\ and 
V2 of the Ba 2+ (2) in the oversized cage vibrating within 
the plane perpendicular to the fourfold axis and along 
the fourfold axis, respectively. It is noted that the Ein- 
stein temperature of about 60 K (=1.2 THz) calculated 
from the atomic displacement parameters is close to 
V\. The broad peak around 2.4 THz can be considered as 
an overlapped spectrum consisting of both the collective 
cage mode V3 and the Ba 2+ (1) mode 1/4 in the smaller 
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FIG. 3: (color online) Temperature dependence of peak fre- 
quencies of BGG for i>i(i = 1 — 6). 



cage. Higher frequencies of 2.98 and 3.30 THz are rather 
close to the collective cage modes v 5 and respectively. 

To extract each phonon contributions <ji^{uj) from the 
total spectra, we subtract the temperature-dependent 
Drude contributions by assuming a linear dispersion. 
Here we adjust the parameters of the linear curve in or- 
der to keep the sum rule that the spectral weight Si of 
Vi is conserved irrespective of temperature; 



ai t i(u>)du = Si = TiNiql/2Mi. 



(1) 



where Mj, Ni, and q\ are the mass, density and charge 
of atoms involved in the V{ phonon mode. Then to the 
phonon spectra from Vi we apply the Lorentz-type con- 
ductivity written by 



oi,i(w) = 



{2l-K)S i u?T i 



(2) 



where Tj is the relaxation rate and wo,i is the resonant fre- 
quency. Here Si is analyzed from Eqfj] using the lowest- 
temperature data having negligibly small Drude contri- 
bution, see Table |TJ The observed ratio of S1/S2 = 1.8 
is close to the theoretical value of 2 taken into account of 
the degrees of freedom for v\ and z/2 modes. Furthermore 
the total contribution (Si + S2) from these Ba(2)-ions' 
modes is about 13.7 x 10 13 , which is larger by factor 4 



than 3.4 x 10 13 calculated with M = M Ba , N = 6/V u 
(V u = 10. 76 3 A 3 is the unit cell volumeQ) and q = 2e for 
Ba(2) ions. Judging from the uncertainty in deriving the 
absolute value of a±, it is not clear at present whether 
this discrepancy is intrinsic or not. To note, however, it 
may be an open issue to study how the charge and mass 
of the guest ion are effectively modified when coupled to 
the encapsulating cage. For 1/5 and v$ modes of cages, 
the small spectral weight as listed in the table may be 
reasonable. 

The peak frequencies o;o,i/27r are plotted as a function 
of temperature in Fig.[3J With decreasing temperature 
down to 6.6 K, the low-lying V\ and ^2 modes soften by 
4.2%, being consistent with the Raman-active modes at 
1.05 THz (35cm" 1 ) with E g symmetry Q, and 1.8% re- 
spectively. On the other hand, other higher-frequency 
modes remain almost constant or harden slightly. In 
general the phonon hardening can be conventionally ex- 
pected due to the lattice contraction at lower tempera- 
tures, so that the magnitude of the softening as observed 
in v\ and ^2 might be underestimated due to an over- 
compensation for the hardening as observed in the cage 
modes v§ and v§. Thus the present studies reproduce 
the softening phenomena in rattling modes of guest ions, 
which have been so far interpreted to be an evident proof 
for an anharmonicity effect in these localized modes. 

In order to shed more light into the anharmonicity ef- 
fect, we focus on the lowest- lying mode v\ to compare 
with model calculations for the optical conductivity. We 
assume the following Hamiltonian describing the motion 
of a single ion in a one-dimensional anharmonic potential 
(1D-AP) including the quadratic and quartic terms for 
the displacement x, 



H 



P 

2M 



1 



kx 2 



1 



Xx 4 



(3) 



where p and M are the momentum and mass (here, sup- 
pose Ba ions). The constants k and A are assumed to 
be positive for a single well. Using a linear response the- 
ory, the optical conductivity &(ui) from the anharmonic 
phonon is expressed as 



a(uj)/ao = iu> /J |(n|x/aro|m)| 

u„„>0 



-i3E„ 



-/3E n 



Z (4) 
x ((w - uj nm + iTo/2)- 1 -{uj + u nm + iTo/2)- 1 ) , 



where \n) and E n are the eigenstate and eigenvalue of 
the Hamiltonian, (3 = l/k B T, Z = J2n e _/3B " and 
k-Vim = (E n — E m )/h. Here a phenomenological relax- 
ation rate Tq is introduced corresponding to the line 
width of the spectral shape. The conductivity is nor- 
malized with do = q 2 Nx 2 ) /%, where Xq is a length scale, 
x o = K/VMk. The calculation indicates that the transi- 
tion probability arises mostly from the matrix elements 
(n + l|x|n) between neighboring levels. Details are dis- 
cussed in Ref. [lfj| . 
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FIG. 4: (color online) (a) Spectral shapes of <ti(lj) from v\ of 
BGG. The filled symbols and solid lines indicate the data and 
calculations respectively. Temperature dependence of (b) the 
peak frequency u>q and (c) the observed line widths compared 
with calculated ones. See the text for the definitions. 



For the numerical calculation of <j\{uj,T) for the v\ 
mode, we fix the parameters: (l/2n)- s /k/M =1.143 THz, 



\/2k x 



= 0.21 (as is Bohr radius), and the 



temperature-independent constant Tq/2tt — 0.067 THz. 
With these parameters the observed spectra are best fit- 
ted to the calculations as shown in Fig. [Ha). Figure. 0Jb) 
shows the temperature dependences of the peak position 
(u ,cai) and the upper (lower) frequencies w u p), ca i giving 
half the peak height of the calculated spectra, together 
with corresponding data Wo. i: andtt^i obtained from 
the Lorentzian fitting. The data can be reproduced with 
our calculations. 

In a usual harmonic approximation with A = 0, E n = 
TiujE{n + 1/2) and all the neighboring level-spacings be- 
come equivalent resulting in the well-known Lorentz-type 
dispersion in the complex conductivity for an harmonic 
oscillator ioj/ (to 2 — w| + iuiTo), as EqO In contrast, the 
anharmonicity due to the quartic term in Eq. Q induces 
unequally spaced energy levels, resulting in n-dependent 
excitations that the larger the n, the wider the spacing. 
At enough high temperatures in comparison to hiOE/ks 
= 55 K (— 1.15 THz), phonons are thermally distributed 
widely from n = to higher levels in accordance to the 
Boltzmann factor, while, at lower temperatures, phonons 
become condensed mostly around n — 0, and can be ef- 
fectively excited to n = 1. The overall spectral shapes, 
given by a superposition of each Lorentzian curves hav- 
ing the n-dependent peak frequencies, become broad at 
higher temperature. Thus the anharmonicity effects on 
the conductivity spectra are expected both in the shift of 
the peak frequency and in the change of the line width. 

Thus the 1D-AP model predicts an asymmetric change 
in the line width with temperature, which can be fea- 
tured in the temperature dependence of 7 tl (i) jCa / = 



l w u(i),ca/ — <^o,cai| measuring the width in the higher 
(lower)-frequency region with respect to the peak fre- 
quency. Figure[4jc) shows stronger temperature depen- 
dence in 7u iC a; than in ji lCa i, indicating the asymmetric 
change of the spectra. The data are well reproduced by 
the mean width defined as half the T ca i — ju.cai + li,cai 
shown by the dashed line. 

Finally, we remark upon two assumptions made in the 
above discussions. First, we have assumed the constant 
line-width (To) in conductivity calculations. This damp- 
ing effect may be attributed, in part at least, to some 
mode-mode couplings to the anharmonic phonons them- 
selves and also to the acoustic cage phonons. Second, 
we have assumed the positive quadratic potential (k > 
in Eq.[3|) meaning the on-center rattling in a single min- 
imum potential. For the case of the off-centering with 
two-level minima with k < 0, the peak frequency and the 
line width are hardly fitted to the data. 

In conclusion, the temperature dependences of the 
peak frequency and the line width in the conductiv- 
ity spectra of BGG are explained by 1D-AP model cal- 
culations without any interactions taken into account. 
Importantly to note, these dependences are determined 
solely by the Bose-factor for optical excitations of anhar- 
monic phonons with the nonequally spaced energy levels. 
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